## ── Attaching packages ─────────────────────────────────────── tidyverse 1.3.1 ──
## ✓ ggplot2 3.3.5     ✓ purrr   0.3.4
## ✓ tibble  3.1.5     ✓ dplyr   1.0.7
## ✓ tidyr   1.1.4     ✓ stringr 1.4.0
## ✓ readr   2.0.2     ✓ forcats 0.5.1
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## x dplyr::filter() masks stats::filter()
## x dplyr::lag()    masks stats::lag()
## Linking to GEOS 3.8.1, GDAL 3.2.1, PROJ 7.2.1
## Reading layer `celdas_marzo_2020' from data source 
##   `/Volumes/GoogleDrive/My Drive/Data/INE Mobility/Supporting Info/shapefiles_ reas_movilidad/Shapefile_celdas_marzo_2020/celdas_marzo_2020.shp' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 3214 features and 6 fields
## Geometry type: MULTIPOLYGON
## Dimension:     XY
## Bounding box:  xmin: -2021660 ymin: 3203377 xmax: 481778.5 ymax: 5433001
## Projected CRS: WGS 84 / Pseudo-Mercator
## Warning in st_centroid.sf(.): st_centroid assumes attributes are constant over
## geometries of x
## Joining, by = "id"
## `summarise()` has grouped output by 'codigo_residencia'. You can override using the `.groups` argument.